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ABSTRACT 

The measurement of the optical depth to free electrons, r r , in the cosmic microwave background 
(CMB) provides an important constraint on reionization, but is degenerate to more complex reion- 
ization models. Small angular-scale CMB measurements of the kinetic Sunyaev-Zeldovich (kSZ) and 
Ostriker-Vishniac (OV) effects have the potential to break this degeneracy. We calculate the OV sig- 
nal for various extended reionization histories described by a simple analytic form. These models are 
parametrized by Az , the duration of the reionization event. For reionization histories with identical 
values of r r , the OV amplitude at I = 3000 (C3000) differs by « 20% between the models with Az = 
0.1 and Az — 3.0. We demonstrate that the removal of the z < 6 component of the OV signal will 
significantly enhance the ability to differentiate between reionization histories, with C3000 varying by 
a factor of ~ 2 between the models with Az = 0.1 and Az = 3.0. If the low-redshift OV and related 
kSZ signal can be adequately subtracted, this would provide an important observational constraint 
on extended reionization models. 

Subject headings: cosmology: theory — cosmology: observations — cosmic microwave background — 
large-scale structure of universe — radiation mechanisms: non-thermal — scatter- 
ing 



1. INTRODUCTION 

The first decade of the 21st century has arguably 
witnessed the advent of precision cosmology, ushered 
in by cosmic microwave background (CMB) experi- 
ments, such as the Wilkinson Micro wave Anisotropy 
Probe (WMAP. iKomatsu et all (|2008f )). and large-scale 
galaxy surveys, such as Sloan Di gital Sky Survey (SDSS, 
lAdelman-McCarthv et all (|2006f) ) and the 2-degree Fie ld 
Galaxy Redshift Survey f2dFORS. IColless et al.l([200lh L 
The major cosmological parameters are now known at 
the percent accuracy level, and ACDM has emerged as 
the cosmology which is most easily reconciled with ob- 
servations. 

The current concordance picture of ACDM cosmology 
holds that after the universe became transparent to the 
CMB at z ~ 1000, the baryonic content of the universe 
was neutral and coupled gravitationally to the dark mat- 
ter component. Some 10 s years after the Big Bang, the 
collapse of massive halos led to the formation of the first 
generation of stars and galaxies. 

These first-generation stars are believed to be ex- 
tremely massive (M > 100M©) and have extremely 
hard spectra. The prodigious amounts of E 13.6 eV 
photons generated by these stars propagated out of 
their host protogalaxies into deep space, gradually re- 
ionizing the inter-galactic med ium (IGM), fully ioniz- 
ing the universe by z ~ 6 (s ee lLoeb fc Barkanal (|2001h 
and iFurlanetto fc Ohl (|2008[ ) for reviews on the sub- 
ject). The reionization process is believed to be highly 
complex: even though considerable theoretical work 
(both analytical and simulation s) have been carried ou t 
to understand th i s epoch (e.g . IFurlanetto et al.1 (|2004f ): 
llliev et all (l2006h ; iTrac fc Cenl ((20071) ), the lack of obser- 
vational evidence has made it difficult to reach consen- 
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sus on fundamental issues such as the temporal e x tent o f 
reionizatio n and when it began (e.g. iLidz et al.l (|2007r ): 
IShin et all (|2008l)), let alone min u tae such as feedback 
proces ses (e.g. ICiardi et al.l (|2000t ): iMesinger &: Dijkstral 
(2008) ) and the det ails of ionized bu bble p ercolation 
(e.g. iMcQuinn et al.l (|2007f k iLee et all (|2008l )) through 
the universe. 

At present, there are two major observational con- 
straints on reionization: (1) The optical depth to free 
electrons up to the surface of last scattering has been 
measured thro ugh CMB polarizatio n data to be r r = 
0.084 ± 0.016 (TDunklev et al J 120081) . corresponding to 
a reionization redshift of z r ps 11 assuming (unrealis- 
tically) that reionization occured instantaneously. (2) 
The detection of the Gunn-Peterson absorption trough 
from neutral hydrogen in the Lyman-a absorption spec- 
tra of high-redshift (z ~ 6) quasars indicat es a rising 
neutr al fraction beginning from this epoch (|Fan et al.l 
2006), hinting at the end of reioni zation process, al- 
though this is by no means certain ([Becker et all 120071 : 
IFurlanetto fc Mesingerll2008l ). 

Over the next decade and beyond, radio experiments 
are being designed to pr obe the 21-cm l ine o f neutral 
hydrogen, e.g. LOFAR (|Zaroubi fc Si lkl 120051) . MWA 
(iMorales fc Hewitt! I2004I ). and SKA (jTerzian fc Laziol 
2006). These 21-cm tomography experiments promise 
full 3D imaging of the neutral IGM, allowing reioniza- 
tion to be studied in unprecedented detail. However, 
these experiments are still in their early stages. 

At time of writing, several high angular resolution 
CMB experiments are ong oing, such as the A tacama Cos- 
mology Telescope (A CT. iKosowskvl (120061) 1 and South 
Pole Telescope (SPT. iRuhl et all (|2004f) ~ These exper- 
iments have resolutions approaching 1', with sensitivi- 
ties of ~ mKy/s. The temperature-temperature angular 
power spectra Ci produced by these experiments will, 
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for the first time, probe small angular scales I > 3000 at 
which the primordial CMB power spectrum is no longer 
the dominant contribution to the anisotropy power. In- 
stead, the secondary anisotropy at these scales will be 
dominated by foreground astrophysics such as unresolved 
radio galaxies, far IR/sub-milimetre dusty galaxies and 
the thermal/kinetic Sunyaev-Zeldovich (tSZ and kSZ, re- 
spectively) and Ostriker-Vishniac (OV) effects. 

In this paper, we study the sensitivity of the OV ef- 
fect towards extended reionization models with the same 
value of r r . While the OV an d kSZ effects have been ex- 
tensively studied in the past (lOstriker fc Vishniac fl986l : 



Vishniad fl 987: Ja ffe fc Kamionkowskil 119981 : iHul l2000l : 
Zhang et al. I~ l2004l) . these studies were made prior to 



the W MAP 5th Data Release (WMAP5, iDunklev et all 
(2008)) constraint on r r (and in the earlier papers, before 
concordance ACDM cosmology was widely accepted). 
iMcQuinn et all (|2005f) made a study on the OV/kSZ sig- 
nal from different reionization models, but this was based 
on preliminary WMAP results, and they did not explore 
the potential of the OV effect to break the degeneracy of 
T r towards different reionization models. 

First, we briefly summarize the derivation of the OV 
effect before discussing our chosen model for parametriz- 
ing reionization histories. We then present our calcula- 
tions and discuss them in the context of upcoming ex- 
periments. 

2. THE OSTRIKER-VISHNIAC EFFECT 

In this section we summarize the derivation of the 
CMB ani sotrop y from the OV effect as elucidated in 
IVishniad (|1987f ) and lJaffe fc Kamionkowskil (1998). The 
fractional temperature change in the direction 9 from the 
scattering of the primordial CMB by bulk motion of free 
electrons is 

at n° r - 

p(9) = = — / n e GT& T {9 ■ v(w9;w)]a(w)dw, (1) 

T Jo 

where n e is the electron number density along the line of 
sight, v(w;w) is the bulk velocity at comoving distance 
w in units of the Hubble distance, at conformal time 77 = 
770 — w (where conformal time is defined by drj = dt/a), 
(Tt is the Thomson scattering cross-section, a(w) is the 
scale factor of the universe at w and we set c = 1. r is the 
optical depth to electron scattering to a given epoch, and 
for a ACDM cosmology it can be expressed as a function 
of redshift thusly: 

n b p c a T a f z x e (z')(l + z') 2 
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where i7{,,f2 m and Q\ are the present-day baryon, mat- 
ter (dark matter + baryon) and dark energy densities 
respectively, in units of the critical density p Cl x e is the 
ionization fraction and m p is the mass of the proton. In 
the second line, we have assumed that x e is constant over 
the period of integration. 
The visibility function is defined as 

dr 

g(w) = n e (w)aTa(w)e T = ——e T . (3) 

dw 

n e from Eq. [1] has been replaced by the mean elec- 
tron number density h e = ribp c x e (z)(l + z) 3 /m, p . g(w) 



represents the probability distribution of first scattering 
from reionized scattering (contrast the analogous visibil- 
ity function in CMB recombination physics, which is the 
probability distribution of last scattering) . The visibility 
function is normalized such that 



g(w)dw = 1 — e 



(4) 



where r r is the total optical depth to scattering by free 
electrons up to the surface of last scattering. 

The fractional temperature change can then be rewrit- 
ten as 

mo 

(5) 



P(0) = ~ / 9{w)9 ■ q(w,w), 



where q(w, w) = [1 + S(w, w)]v(w, w) and <5(w, w) is the 
fractional density perturbation. 

Using the linear theory approximations for <5(w, w) and 
v(w, w), and projecting the temperature perturbation 
onto the sky using the Limber approximation, we get 
the Ostriker-Vishniac angular power spectrum in terms 
of multipole moments 



a = 
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where w = 770 — f], D is the growth factor (the dotted 
variable denotes a derivative with respect to time) and 
the function S(l/w) (sometimes known as the Vishniac 
power spectrum) is usually written in terms of the co- 
moving wavenumber fc: 



S(k) = k dy dx P(ky)P(ky/l + y 2 
Jo J -1 

y {l-x 2 ){l-2xyf 

(1 + y 2 - 2xy) 2 



2xy) 



(7) 



where P(k) is the standard linear theory matter power 
spectr um using the transfer function from Bardcc n et al.l 
(1986), normalized to the latest WMAP5 values. 

We now make a brief note on nomenclature regard- 
ing the related kSZ and OV effects, which have at 
times been used interchangeably by authors. The tem- 
perature perturbation equatio n (Eq. []} was first de- 
rived by iZeldovich fc Sunvaevi (|1969l ). and then stud- 
ied in the context of t he peculiar motion o f ind i- 
vidual galaxy clusters in ISunvaev fc Zeldovichl (| 19801) . 
lOstriker fc Vishniad (|1986h and lVishniad (|1987l ) were the 
first to study the effect in the context of large scale 
structure, specifically with linear perturbation theory . 
We therefore refer to the OV effect as the temperature 
anisotropy caused by velocity and density fields in the 
linear regime, while the kSZ is its non-linear counterpart 
(which is the effect that accounts for the temperature 
perturbation from the peculiar motion of galaxy clus- 
ters). 

3. REIONIZATION HISTORIES 

One of the most important outstanding questions on 
reionization is: when did it happen? 

The naive way of answering this is by measuring the 
total optical depth r r to free electrons from us to the 
surface of last scattering. One can then invert Eq. [5] 
to find the epoch of reionization z r . This estimate of 
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Fig. 1. — Ionization fraction x e (left) and visibility function g (right) as a function of redshift, for 3 reionization histories constrained 
by t> = 0.09. The discontinuity at 6.0 < z < 6.5 is caused by our linear interpolation scheme which ensures that x e (z = 6.0) = 1 (see 
discussion in Sec.[3]l. 

Prior analyses of the kSZ or OV effects have tended to 
assume instantaneous reionization, in which the visibility 
function g(z) has been calculated assuming that x e = 1 
at z ^ z r and x e = prior to that. 

In order to describe an extended reionization event, 
we choose a simple parametrization for the mean neutral 
fraction as a function of redshift, ofte n used to compare 
resul t s of different reioniz ation models (fZaldarriaa a et all 
2004; I Alvarez et aHl200l 




Fig. 2. — Reionization histories parametrized by z r and Az (as 
described by Eq. [g}, solved for r r = 0.100 (solid line), r r = 0.085 
(dotted line) and r r = 0.070 (dashed line). 

z r assumes an instantaneous reionization event which is 
physically unrealistic, while r r by itself is degenerate to 
more complicated reionization histories. 

Until recently, measurements of r r have had consider- 
ably worse precision, and authors studying reionization 
have tended to regard r r as a tunable free parameter in 
their reionization models. The WMAP5 measurement 1 of 
T r = 0.084 ± 0.016 (corresponding to an instantaneous 
reionization redshift z r — 10.8 ± 1.4) has been precise 
enough to change this paradigm, and in this paper we 
regard r r as a fixed value, albeit one that still has non- 
negligible uncertainties. However, one can expect sub- 
sequent WMAP data releases as well as the upcoming 
Planck space mission (|The Planck Collaboration] 2006) 
to further tighten the constraint on r r to < 10% preci- 
sion. 



xhi(z) 



1 



1 + exp[— (z — z r )/ Az] ' 



(8) 



where z r in this case is defined as the redshift at which 
%Hi = 0.5 and Az parametrizes the duration of the reion- 
ization event. The mean ionized fraction is then trivially 

X e = 1 - XHI- 

Up to this point, we have assumed that the universe is 
composed solely of hydrogen, and ignored helium reion- 
ization. It is often assumed that the He I in the universe 
becomes singly ionized at roughly the same time as H I 
as their (first) ionization energies are of the same mag- 
nitude (24.6 eV for He I compared with 13.6 eV for H I). 
To make this assumption, we modify the expression for 
the number density of electrons used in Eqs.[T][3] 



1 -Y 



Y 

4m r , 



Q, b p c X e (z){\ + zf 



(9) 



1 http://lambda.gsfc.nasa.gov 



where Y is the helium mass fraction of the universe (in 
our calculations we use the value Y = 0.24). The two 
terms in the square parentheses correspond to the free 
electrons contributed by the ionization of hydrogen and 
the first ionization of helium, respectively. 

He II reionization, on the other hand, is a low-redshift 
(z ~ 3) event which is much more amenable to observa- 
tions (jShull et al.H2"00llFurlanetto fc Ohll2"008l ) and thus 
beyond the scope of this paper. 

While the functional form of Eq. [8] is believed to de- 
scribe well the evolution of the neutral fraction during 
most of the reionization process, it becomes less accurate 
near the end of reionization as xri 0. For reionization 
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Fig. 3. — Angular power spectra from the Ostriker-Vishniac effect, plotted logarithmically (left) and linearly (right). The solid, dotted 
and dashed lines correspond respectively to reionization histories with Az = 0.1, 1.0 and 3.0. 
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Fig. 4. — Angular power spectra from the Ostriker-Vishniac effect calculated for z > 6, plotted logarithmically (left) and linearly (right). 
The solid, dotted and dashed lines correspond respectively to reionization histories with Az = 0.1, 1.0 and 3.0. Note that the scale on the 
abscissa differs from that in Fig. \3\ 



histories consistent with r r w 0.09 and z r w 10 the neu- 
tral fraction doesn't drop quickly enough to Xhi ~ 10~ 3 
by z < 6.5, as observed in high-redshift quasar obser- 
vations. We thus modify our reionization model as fol- 
lows: the neutral fraction is parametrized by Eq.[5]in the 
regime z > 6.5, and then interpolated linearly to xhi = 
at z — 6.0. The left-hand plot of Fig. [1] illustrates our 
chosen form of x e . 

For a fixed T r , this gives a family of reionization mod- 
els parametrized by z r or Az . In Fig. [2] we compute Az 
and z r for reionization histories constrained by 3 different 
values of r r consistent with WMAP5 (hereafter we use 
Az to label different reionization histories at fixed r r ). 
We see that any given value of r r is degenerate with a 
wide range of reionization histories ranging from instan- 
taneous (Az — > 0) to highly extended (Az ~ 4—5). Note 



that smaller values of r r tend to limit the possible range 
of Az due to our constraint that reionization completes 
by z ~ 6, as discussed above. 

Fig. Q] shows our chosen form of the ionization frac- 
tion x e and the corresponding visibility function g as a 
function of redshift, calculated for several reionization 
histories which share the same value of r r . 

4. ANGULAR POWER SPECTRUM CALCULATIONS 

Eqs.[6]and[7]are integrated numerically to calculate the 
OV angular power spectrum. In ou r calculations we use 
the WMAP5 recommended values (jDunklev et al.ll2008| ) 
for the cosmological parameters, h — 0.71 (where Hq = 
100ft. km s^Mpc" 1 ), n m = 0.28, fl b h 2 = 0.0227, n = 
0.96, A 2 R (k = 0.002 Mpc- 1 ) = 2.41 x 10~ 9 . 

However, dark energy drops to less than 10% of the 
total mass-energy of the universe above z ~ 1.3, while 
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the main epoch of interest at which reionization affects 
the OV signal is z > 6. We thus make the simplifying 
assumption of a flat Einstein-de Sitter universe (J7a = 0, 
= 1) for computing the OV angular power spectrum. 
This allows us to make the substitutions a — rj 2 frfa and 
DD/Dq = 277/770 in Eq. [6] The integrated optical depth 
r(z), however, is highly sensitive to cosmology and we 
assume ACDM cosmology when calculating this for the 
visibility function (c.f. Eq. [3]). 

While comparing the results of ou r calculations with 
those in lJaffe fc Kamionkowskil (| 19981 ). we found discrep- 
ancies that were eventually traced to an error in their 
code for calculating S(k) (c.f. Eq. [7])- This error caused 
an overall underestimate (by ~ 10 — 20%) in the OV 
power of the results presented in their paper. The shape 
of their l 2 Ci power spectra should also be slightly more 
peaked, and the peaks have tended to move to lower I 
after the mistake was corrected. Hence, the analytic ap- 
proximation for S(k) fitted to their results (Eq. 41 in 
their paper) is incorrect. Nevertheless, the equations de- 
rived in their paper are mostly correct, except for their 
Eq. 38, which should be identical to the second line of 
our Eq. [2 

The resultant power spectra for several reionization 
histories consistent with r r = 0.09 are shown in Fig. [3l 
The power spectra peak at I ~ 1000 while the variance 
of the temperature distribution 

{^P)) = hj lCldl < 10 > 

corresponds to a root-mean-squared (r.m.s.) value AT ~ 
2.4 — 2.7 /iK. The amplitude of the spectra tends to 
decline with increasing Az, as this has the effect of 
spreading g(z) to higher redshifts (c.f. Fig. [T|) while the 
(aDD /wDq) 2 parts of the integrand are decreasing func- 
tions with redshift. Hence, a more extended reionization 
event samples the low-redshift part of the integral less 
than a rapid event, which decreases the power in the 
signal. 

At I — 3000, where the primary anisotropy drops be- 
low the level of the secondary anisotropies, C3000 differs 
by about 20% between the cases of near-instantaneous 
(Az = 0.1) and highly extended (Az = 3.0) reionization. 

However, for the majority of the universe's lifetime we 
know the ionization state: x e = 1 from z < 6 onwards 
irrespective of the reionization model, so if we remove the 
low-redshift part of the integral we get rid of the part of 
the signal that is independent of the reionization model. 

In Fig. [4] we plot the contribution of different redshifts 
in the calculation of C3000 (i-e. the integrand of Eq. [6] 
plotted against redshift for I = 3000). We see that a sig- 
nificant portion of the power is contributed from z < 6. 
Removing this contribution would increase the contrast 
between the Ci 's computed for different reionization his- 
tories. 

This is duly shown in Fig. [3l where we have removed 
the z < 6 portion from the integral in Eq. [6l The re- 
maining power in the anisotropies peak at I ~ 1500, 
but now the r.m.s. temperature perturbation varies from 
AT = 1.5 /iK for Az = 0.1 to AT = 1.0 jiiK for Az = 3.0. 
C3000 n ow varies by a factor of ~ 2 between the near- 
instantaneous reionization model and the most extended 
one. 
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Fig. 5. — OV anisotropy at I = 3000, plotted as a function of the 
contributing redshift for Az = 0.1 (solid), Az = 1.0 (dotted) and 
Az = 3.0 (dashed). 

5. DISCUSSION / CONCLUSION 

From our calculations, we have shown that by altering 
the duration of the reionization event that goes into the 
visibility function of the Ostriker-Vishniac power spec- 
trum, significant changes are made to the resulting C;'s. 
However, the OV effect is certainly not the only sec- 
ondary anisotropy to exist at I > 10 3 . 

Assuming that instrumental noise issues have been 
dealt with, the main challenges facing the detection 
of the kSZ/OV signal at I > 2500 are the other sec- 
ondary anisotropies. The galactic foregrounds and far- 
IR/sub-mm extragalactic sources can be removed by 
cross-correlating with multi-wavelength data, while the 
thermal SZ signal can be removed due to its spectral 
dependency near the CMB peak temperature. 

In this paper we have only considered the linear OV 
effect, and not the full combination of OV (which dom- 
inates at high-z) and non-linear kSZ (which dominates 
at low-z). However, the calculation of kSZ is carried 
out through an integral similar in form to Eq. [6l with 
the main differences being the calculation of the inte- 
gral S(k) (Eq. [7]). In practice, these involve using some 
form of non-linear po wer spectrum instead of th e linear 
P(k) we used in Eq.[7l pIull2000l : lMa fc Fry|[2002]) . These 
calculations are still weighted by the electron visibility 
function g(z) in the same manner as the OV, but the 
fact that the non-linear terms are relatively small com- 
pared to the linear terms at z > 6 mean that the kSZ 
is l ess sensitive to the reionization histo ry than the OV; 
see lHernandez-Monteaeudo fc Hoi (|2009f ) for a recent pa- 
per on this. One additional source of OV/kSZ power 
is from patchy reionization, i.e. the contribution from 
large ionized H II bubbles which form during reioniza- 
tion ([McQuin n et al.ll2005f ). However, this contribution 
should similarly be weighted by the visibility function 
and would not qualitatively affect the dependency of the 
kSZ/OV signal on the reionization history as discussed 
in this paper. 

In the most optimal case, the noise sensitivity of ACT 
would allow for Ci measurements at a precision of sev - 
eral percent at I of several thousand (|Kosows 
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Once the systematic errors from the subtraction of fore- 
grounds, extragalactic point-sources and tSZ are taken 
into account, it would be challenging to measure the full 
kSZ/OV signal to sufficient accuracy to distinguish be- 
tween different reionization histories, for which the val- 
ues of C3000 var y by about 20% between the cases of 
instantaneous reionization and highly extended (Az = 
3) reionization. 

The key to distinguishing between reionization histo- 
ries lies in removing the low-redshift (z < 6) component 
from the kSZ/OV signal (specifically the kSZ contribu- 
tion from non-linear structures, which dominate at low 
rcdshift). This would greatly increase the contrast be- 
tween the different reionization histories, to a factor of 
~ 2 between the cases of Az — 0.1 and Az = 3.0 as we 
h ave seen above. 

iHo et all (|2009h have developed a method of measur- 
ing the kSZ signal at low redshift. They build a template 
of the momentum field by reconstructing the momen- 
tum field from galaxy redshift surveys, and then cross- 
correlating this with arcminute-scale CMB data. For a 
cross-correlation between a 4000 sq. deg. ACT exper- 
iment with galaxy survey data from the 3rd phase of 
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the SDSS (SDSS3), IHo et all ((2001 estimate a kSZ de- 
tection with a signal-to-noise (S/N) of ~ 15. However, 
this method can only measure the kSZ contribution from 
relatively low redshifts (z < 1) that can be adequately 
covered by the galaxy redshift surveys of the foreseeable 
future. 

One possibility would be to construct a kSZ tem- 
plate from the Lya forest in the spectra of high-redshift 
quasars, which coincidentally covers the z < 6 redshift 
region for which we want to subtract the kSZ. This would 
be complicated by the difficulty of extracting the under- 
lying density field from Lya data, but even a relatively 
low S/N estimate of the low-redshift kSZ should allow us 
differentiate between instantaneous and highly extended 
reionization. Considering the current paucity of observa- 
tional data on reionization, even a rough estimate of Az 
would provide valuable constraints on reionization. 
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